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Abstract 

The thermodynamic properties of a SU (3) gauge theory without quarks are calculated 
using a string formulation for 1.2 T c <T< 3T C . The results are in good agreement with the 
lattice data. We also comment on SU(N) gauge theories. 
PACS numbers: 12.38.Lg, 12.90.+b 



1 Introduction 

High energy nucleus-nucleus collisions provide the means of creating nuclear matter in conditions 
of extreme temperature and density. In particular, the system undergoes a transition from a state 
of nucleons containing bound quarks and gluons to a state of deconfined quarks and gluons. This 
state was originally given the name Quark Gluon Plasma. However, the results at RHIC indicate 
that instead of behaving like a weakly coupled gas of free quarks and gluons, the matter created 
in heavy ion collisions behaves like a strongly coupled liquid^ Thus, there is a need for new 
approaches to strongly coupled gauge theories. 

Until recently, the lattice formulation was a unique theoretical tool to deal with strongly 
coupled gauge theories. The subject has taken an interesting turn with Maldacena duality [2]. 
One of the implications is that it resumed interest in finding a string description of strong inter- 
actions. Although the original proposal was for conformal theories, various modifications have 
been found that produce gauge/string duals with a mass gap, confinement, and supersymmetry 
breaking [3]. 

In this paper we address some issues of thermodynamics of SU (3) pure gauge theory in a dual 
formulation. Clearly, finding the dual from first principles of string theory is beyond of our ability. 
Instead, we attempt the inverse problem and use our knowledge of some phenomenologically 
successful five-dimensional models of AdS/QCD. 
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Before proceeding to the detailed analysis, let us set the basic framework. We consider the 
following ansatz for the 10-dimensional background geometry which turns out to be applicable 
for the temperature range 1.2 T c < T < 3T c |j 

ds 2 = ^H(fdt 2 + dx 2 + f- 1 dz 2 )+H- 1 g ab dio a cLu b , f = 1 - , H = e^ 2 ' . (1) 

Here z T = 1/nT. It is a deformed product of the Euclidean AdSs black hole and a 5-dimensional 
sphere (compact space) whose coordinates are u) a . The deformation is due to a z-dependent factor 
H. Such a deformation is crucial for breaking conformal invariance of the original supergravity 
solution and introducing A QCD . 

Apart from the language of 10-dimensional string theory, there is a more phenomenological 
way to attack QCD. This approach called AdS/QCD deals with a five-dimensional effective 
description and tries to fit it to QCD as much as possible. For our model, its AdS/QCD cousin 
can be obtained by discarding the compact space in (P). 

At T = 0, then what we get is the slightly deformed AdS5 metric. Such a deformation is 
notable. The point is that in this background linearized Yang-Mills equations are effectively 
reduced to a Laguerre differential equation. As a result, the spectrum turns out to be like that 
of the linear Regge models [U [5] . This fact allows one to fix the value of c from the p meson 
trajectory. It is of order [5] 

c«0.9GeV 2 . (2) 

We will assume that the value of c is universal and is therefore valid for the world without quarks 
too. In addition, this AdS/QCD model provides the phenomenologically acceptable heavy quark 
potentials as well as the value of the gluon condensate [U [7]. 

At finite T, the model provides the spatial string tension of pure gauge theory [8]. The 
agreement with the lattice data is very good for temperatures lower than 2.5-3 T c . Due to this 
reason we set the upper bound on T in (JT() . Moreover, the model describes in a qualitative way 
a heavy quark-antiquark pair and the expectation value of the Polyakov loop [9] . 

Thus, there are reasons to believe that the model ([I]) is a good approximation for a string 
dual to a pure gauge theory. 

2 Thermodynamics 
2.1 The Entropy Density 

One of the bedrocks of gauge/string (gravity) duality is a conjecture that the entropy of gauge 
theories is equal to the Bekenstein-Hawking entropy of their string (gravity) duals [3]. As known, 
the Bekenstein-Hawking entropy is proportional to an (8-dimensional) area of the horizon. We 
can now ask whether the five-dimensional framework (AdS/QCD) is an adequate approximation 
at this point. In general, the answer is no. There is a contribution from the compact space that 
might be relevant. 

2 The lower limit is chosen to keep the system out of the critical regime. As we will discuss below, the upper 
limit is determined by consistency rather than perturbation theory. 
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The metric (Q"]) has the horizon at z = z T . Therefore, the temperature dependence of the 
entropy density iqf| 

s(T) = s T 3 exp{-^}, (3) 
where so is a factor independent of temperature. In this formula T c is given b)@ 

T c = -V~c (4) 

7T 

It follows from ([3]) that the entropy density can be represented as a series in powers of ^ 
with the leading T 3 term 

oo „ 2 

s(T)=s T 3 ^2a n r n , r=^, (5) 

n=0 

(—) n 

where a n = 

For future use, we define the truncated model by keeping the two leading terms in ©. We 
have 

Str (T)= S0 T 3 (l-ir). (6) 



2.2 The Pressure 
2.2.1 Fuzzy Bags 

Recently, it has been suggested by Pisarski that for the temperature range T max < T < T pcrt the 
pressure in QCD with quarks is given by a series in powers of ^ times the ideal T 4 term |10j . 
Explicitly, 

Pqcd(T) ~ f pcrt T — i3f uzzy T — B MiT + . . . . (7) 

It was called a fuzzy bag model for the pressure. So, B MIT stands for the MIT bag constant. 
T max is close to a critical temperature T c (or some approximate "T" for a crossover). A small 
difference between T c and T max may vary with the model. T pert is set by perturbation theory such 
that it is applicable only for temperatures higher than T pcrt . 

For pure glue, Pisarski argued, based on lattice simulations of [11 j . that ([7]) reduces to 

p(T)^f pcit (T 4 -T?T 2 ) . (8) 

This means that B fuzzy = f pcit T 2 and B MIT is much smaller than the first two terms. So, the 
pressure is a sum of just two pieces. Note that an important consequence of ([8]) is that the 
pressure (nearly) vanishes at T = T C . 

3 We take a constant dilaton. 

4 In the following section we will see that T c can be thought as a critical temperature. 
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2.2.2 String Dual 



Given the entropy density as a function of T, in the homogeneous case one can find the temper- 
ature dependence of the pressure by integrating 4& = sH From ©, we get 



i 1 00 

p(T) = -s T 4 (l- T --T 2 \nT-bT 2 + Y. 6 « rn ) > 



(9) 



n=3 



where 6 is an integration constant and b n = 



The final topic to be considered here is whether the proposal of Pisarski is reasonable in the 
model under consideration. The two leading terms in ([9]) look similar to those of (jHJ). So, we 
find that the critical temperature is given by T c . A simple estimate then givea^] 



T r w 300 MeV . 



(10) 



In SU(3) pure gauge theory the critical temperature is of order 270 MeV. So, the agreement is 
not bad at this point. 

Let us now use p(T c ) = to determine the integration constant. As a result, we have 



6 = ^6„« 0.039. (11) 

n=3 

The value of b is indeed small compared to the coefficients in front of the two leading terms. 
Thus, the agreement is very satisfactory at this point. 

To complete the picture, we present the results of numerical calculations. We split the series 
([9]) into two pieces, the first containing the two leading terms, and the second presenting the 
rest. Then we defin^] 



Pi (T) = 1-t. 



P2(T) 



1 x 

-t 2 lnr - br 2 + ^ b r , 



(12) 



n=3 



For the sake of simplicity, we have omitted the overall factor jSqT 4 ". The values of p\ and p2 
can be read off of Fig.l. We see that at T ~ 1.2 T c the value of p2 is one order of magnitude 




1.25 1.5 1.75 



2.25 2.5 2.75 3 



TIT, 



Figure 1: Values of p\ and p2 versus the ratio ^r. 



5 In what follows, we consider the homogeneous case. 
6 We use ((2)1 for this estimate. 

7 Note that the truncated model ([5| can be derived from pi times ^so T 4 . 
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smaller than that of p\. Above 1.2 T c the value of p\ increases, while P2 decreases and becomes 
negligible for T > 2T C . Thus, pi(T) provides a reliable approximation whose error is less then 
10% for the pressure. 

In sum, the truncated model which is equivalent to the proposal of Pisarski is valid with 
accuracy better than 10%. 

2.3 The Speed of Sound 

Having derived the entropy density, we can easily obtain the speed of sound. For the model of 
interest, we have 

For completeness, we also present the result obtained for the truncated model ©. In this case 
(p~3|) is replaced by 

^^(l-irKl-ir)-'. d4) 

Note that C s is independent of so. Thus, we do not have any free fitting parameter at this point. 

We close the discussion of the speed of sound by comparing the results with those of lattice 
simulations H The curves are shown in Fig. 2. 

c 2 s 



0.32 




Figure 2: The square of the speed of sound versus jr. The upper and lower dashed curves correspond 
to (TT5|) and (TTJ} , respectively. The solid curve represents the result of the extrapolation to the continuum 
limit for lattice simulations [IT] . The solid horizontal line is the usual AdS/CFT result with the value |. 



From the Figure we can see that the model of interest is in very good agreement with the 
lattice for T > 1-7T C , while near 1.2T C the discrepancy is of order 15 percent. The agreement 
between the truncated model and the lattice is spectacular. The maximum discrepancy occurred 
at T = 1.2 T c is of order 6 percent. 

2.4 The Gluon Condensate at Finite Temperature 

We will next describe the gluon condensate at finite temperature^ It is obtained from the trace 
anomaly of the energy-momentum tensor [13]. We have 

8 The recent data of [12] have large error bars. So, it is impossible to say how precisely the results fit. 
9 Although the literature on the gluon condensate is very vast, to our knowledge, there are no reliable results 
for the temperature range 1.2 T c < T < 3T C except those coming from lattice simulations. 
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G 2 {T) = G 2 + 4p-Ts, (15) 

where G 2 is the condensate at zero temperature. 

Unlike the speed of sound, the condensate depends on the parameter sq. There are two 
different ways to fix its value which fortunately yield very similar results. The first is to fit the 
interaction measure (e — 3p)/T 4 as it follows from ([5|) and ([9]) to the lattice data of [11] at some 
normalization point T n . As result, we get 

s = 6.8 ±0.3. (16) 

At first glance it may seem curious that the result is almost independent of the normalization 
point. As we will see in a moment, this is indeed the case. 

The second is to match the coefficient in front of the T 4 term in ([9]) with that of the bag 
model [14] . For SU(N) (pure) gauge theory, the latter is simply jV 4 ^" 1 7r 2 . At N = 3, we find 

so = §^7.0 (17) 

that is really the same as (|16p . 

Having determined the value of so, we can now write down the expression for the condensate. 
Combining ([15]) . and © and ©, we get 

G 2 (T) = -s T 4 Qr + - A r 2 In r + gr 2 + f> n - b n )r^j , (18) 

where g = | + b— ^. Note that the condensate at zero temperature G 2 = kT 4 has been included 
in the r 2 term. For the background geometry ([I]), the estimate of [7] gives k ~ 1.20. Interestingly, 
the value of g turns out to be small. For so = 6.8 it is of order —0.01. 

In Fig. 3 we have plotted the gluon condensate in units of T 4 as a function of the ratio jr. 

Gift 

T/T c 
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Figure 3: The gluon condensate in units of T 4 versus Sr. Here sq = 6.8. 



We conclude the discussion with a couple of comments: 
(i) The expression for the gluon condensate is cumbersome and difficult of any practical use. 
We should therefore seek a simpler (nearly equal) expression. What we already know is that 
the truncated model is a good approximation for the pressure and the speed of sound. So, it is 
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reasonable to use this option. We can check it by the same arguments that we used in the case 
of the pressure. To this end, we split the series (fT8j) into two pieces and define 

1 1 00 

91 (T) = --r , g 2 (T) = --t 2 lnr - gr 2 + J> n - a n )r n . (19) 

n=3 

For simplicity, we have omitted the overall factor sqT 4 . 

The values of g\ and g 2 can be read off of Fig.4. We see that the value of g 2 is approximately 

9 2 
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Figure 4: Values of g\ and 32 versus the ratio 



15% of g\. Thus, in the temperature range under consideration we may approximate the infinite 
series (fTHj) by g\ . Finally, the gluon condensate takes the form predicted by the truncated model 

G 2 (T)^-^T 2 T 2 . (20) 
(ii) Using ([5]) and Q, one can easily find the expression for the interaction measure. It is 
e — 3p 



rp4 



(11 1 \ 

_ r + _ T 2 ]nT + ( 6 + _) r 2 + J2(a n - b n )r n ) . (21) 



n=3 

The truncated model provides a simpler expression of the measure 

^ = fr, (22) 

as expected. 

In Fig. 5 we have plotted the interaction measure as a function of the ratio jr. As can be 
seen from the Figure, the agreement with the lattice data is very satisfactory. An important 
observation is that varying sq over the range (|16h has a little effect. 



3 Concluding Comments 

(i) The model we have proposed predicts the entropy density as a series in It differs from 
the proposal of Pisarski [10] by having a term In T in the pressure. However, in the pure glue 
case the InT term turns out to be subdominant. 

(ii) Interestingly enough, the spatial string tension calculated within the AdS/QCD cousin of (PQ) 
is given by [8] 
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Figure 5: The interaction measure (e — 3p)/T 4 versus Sr. The red and green curves correspond to (|2ip 
and (T21Zj) . respectively. The lattice data of [TT] lie between the two dashed lines. Here sq = 6.8. 



e- 1 



a s = a , (23) 

T 

where a is the string tension at zero temperature. Then, from (|23p . we see that the spatial tension 
can be written as a series in powers of ^ times T 2 . Note that unlike the cases considered in 
section 2 the first two terms of the series do not provide a reasonable approximation, 
(hi) Can one think of the model (pQ) as a string dual to a SU(N) pure gauge theory? We will be 
exploring the consequences of assuming that the pressure vanishes at T = T c 

This assumption leads to the same expression for the pressure as Q with b defined by (fTT|) . 
The overall constant so is fixed from the T 4 term. Fitting the bag model, we have 

so = ^ (iV 2 - 1) • (24) 

Clearly, the analysis of section 2 is not sensitive to N. So, the conclusion we draw is that the 
truncated model which is equivalent to dH) is valid with accuracy better than 10%. 

Moreover, we can obtain a formula for the pressure normalized by the leading term po = 
\s T 4 . It is 

^ oo 

-P(r)= (l- T --T 2 \nT-bT 2 + Yb ri T n ). (25) 

Thus our model predicts that the ratio is a function of ^-0 It does not explicitly depend on 
N . At this point it is worth noting that in addition to iV = 3 the prediction is also supported 
by lattice simulations for iV = 4 and N = 8 [15j . 

(iv) We can gain some understanding of the N dependence of a parameter g = Here a' is 
the usual string parameter coming from the Nambu-Goto action. 
The lattice data are well fitted by [16] 



-| = 0. 5 % + ^, (26) 



10 The parameter c is now dependent of N, so is not applicable. 

"1 1 

Strictly speaking, it is a function of =§. 
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where a is the string tension at zero temperature. For the AdS/QCD cousin of ([I]) it is given by 




. = 5-0. (27) 



Combining ([26]) . and (J2|) and ([271) . we learn 

4 / 0.453 V 2 , . 

A simple algebra shows that q is a slowly varying function of N. It takes values between 0.93 at 
N = 2 and 1.32 at N = oo. 

For N = 3, q is approximately equal to 1.12. It is interesting to compare this value with the 
estimate of [6]. The latter was made by using the Cornell potential. The result is g ~ 0.94. The 
estimates are relatively close. This might be a hint that q is also a slowly varying function of a 
number of quarks. 
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